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Dynamical Jahn-Teller effect has escaped so far direct observation in metallic systems. It is 
particularly believed to be quenched also in correlated conductors with orbitally degenerate sites 
such as cubic fullerides. Here the Gutzwiller approach is extended to treat electron correlation over 
metals with Jahn-Teller active sites and applied to the investigation of the ground state of KaCeo- It 
is shown that dynamical Jahn-Teller instability fully develops in this material when the interelectron 
repulsion Lf on Ceo sites exceeds some critical value. The latter is found to be lower than the current 
estimates of Lf, meaning that dynamical Jahn-Teller effect takes place in all cubic fullerides. This 
leads to strong splitting of LUMO orbitals on Ceo sites and calls for reconsideration of the role of 
orbital degeneracy in the Mott-Hubbard transition in fullerides. 

PACS numbers: 71.70.Ej, 71.20.Tx 71.27.+a 


I. INTRODUCTION 

Dynamical Jahn-Teller effect (JTE) is an ubiquitous 
phenomenon in molecules and isolated impurity centers 
with orbital degeneracy.ll^^ Its presence in Jahn-Teller 
crystals is encountered less often, where cooperative or¬ 
dering of static Jahn-Teller distortion is the most prob¬ 
able scenario.l^ Dynamical JTE has been advocated as a 
reason for the lack of orbital ordering in some insulating 
materials such as LiNi02|^Ba3CuSb209,l^and EeSc2S4.^ 
It was also assumed to take place in insulating fullerides 
AiCeo with A = K, Rb, Cspand Li 3 (NH 3 ) 6 C 6 o!^ Re- 
cently, ab initio calculations have shown that dynami¬ 
cal JTE is the reason for the lack of orbital ordering in 
CssCeo fullerides, which explains their conventional an¬ 
tiferromagnetic ordering.^ As for metallic systems, no 
direct evidence for development of dynamical JT insta¬ 
bility in their ground state has been obtained so far. 
Materials where such instability is likely to be realized 
are metallic cubic fullerides A3 Ceo for the reason that 
these correlated conductors are close to Mott-Hubbard 
insulator^^^^hJI for which the existence of dynamical JTE 
was already proved 

In A3 Ceo, the conduction band originating from 
the triply degenerate tiu lowest unoccupied molecu¬ 
lar orbitals (LUMO) on fullerene sites strongly cou¬ 
ples to the i ntram olecular JT active fivefold degener¬ 
ate hg modesDespite the JT coupling, the sym¬ 
metry lowe ring has not been observed in x-ray diffrac¬ 
tion datgP^lh J I ^ I ^^ I implying that the JT effect is either 
quenched by the formation of the band or dynamical. 
An adequate description of the ground state in metallic 
^3 Ceo requires a concomitant treatment of the JT effect 
and the electron correlation. One of the simplest meth¬ 
ods to treat the electron correlation is by v ariatio nal ap¬ 
proach with the Gutzwiller’s wave function.^^^^^ Despite 
the simplicity, Gutzwiller’s approach allows to take into 
account the main contribution to the correlation energy. 
Concerning the ground state of metallic phase, the de¬ 
scription by this method is comparable in accuracy to 
dynamical mean-field theory.!^ Moreover, it has been ex¬ 


tended to treat various situations, for instance, the multi¬ 
band systems . 1 ^ However an adequate approach suitable 
for degenerate conductors with JT effect on sites is still 
lacking. 

In this work, we propose a method to treat elec¬ 
tron correlation in metallic JT systems based on a self- 
consistent multiband Gutzwiller ansatz and apply it to 
metallic fullerides. We find that a dynamical JT instabil¬ 
ity takes place in A3 Ceo already at intermediate strength 
of electron correlation, leading to large amplitudes of JT 
distortions on the fullerene sites and to the removal of 
the degeneracy of three LUMO levels. This means, in 
particular, that the electron correlation in A3 Ceo devel¬ 
ops not in a threefold degenerate LUMO band as thought 
befor^^but in three split subbands. The immediate im¬ 
plication is that the degeneracy of the LUMO band as 
a reason for the h igh c ritical value Lf/w for the Mott- 
Hubbard transitiorP^l^ {w is the bandwidth of the de¬ 
generate LUMO band) should be reconsidered for ful¬ 
lerides. 


II. ELECTRONIC AND VIBRONIC MODEL 
FOR THE LUMO BAND IN METALLIC 
FULLERIDES 


The model Hamiltonian of A3 Ceo consists of trans¬ 
fer Hi, Jahn-Teller L^jt, and on-site bielectronic 
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FIG. 1. (color online) Orientation of Ceo with respect to 
tetragonal axes of fee lattice. 
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where m is a site, Am is a position relative to m, 
A, A' = z are the components of the tiu LUMO (Fig. 
[^, cr is the spin projection, 7 = is the com¬ 

ponent of the hg vibrational mode (7 = 1,4, 5, 2,3 in 

Ref. 1271 respectively), cJ^Acr(^mAcr) is the creation (anni¬ 
hilation) operator of an electron in orbital Act at site m, 

^mAcr = cJnAcr^mAcr^ ^nd q^'y and pjn 7 are the dimension¬ 
less normal coordinate and its conjugate momentum,^ 
respectively. is the Clebsch-Gordan coefficient 

and g are the frequency and the dimensionless vibronic 
coupling constant for the effective hg mode, is the 
transfer parameter, U\\ and U±_ = U\\— 2Jh are the intra 
and interorbital Coulomb repulsion on the fullerene site, 
respectively, and Jr is the Hand’s rule coupling. 

The tight-binding Hamiltonian © has been 
parametrized on the basis of density functional theory 
(DFT) band structure calculation of K 3 C 60 and includes 
nearest neighbor and next-nearest-neighbor electron 
transfer (see Appendix [A| for details). Although nearest 
neighbor tight-binding models were intensively us ed in 
the past to describe the LUMO bands of fnlleridesj ^^FQ I 
the inclusion of next-nearest-neighbor electron trans¬ 
fer is necessary for realistic description of the band 





FIG. 2. (color online) The coordinate systems used to de¬ 
scribe JT effect on a fullerene site. x,y,z correspond to the 
orthorhombic LUMO orbitals (Fig. and 1, 2, 3 to the adia¬ 
batic orbitals. x',y',z' and x",y",z are intermediate coor¬ 
dinate system^^ appearing during the Euler rotation of the 
orbitals from x,y,z to 1,2,3, Eq. 


dispersion. F^Fd xhe JT effect in fullerene anions in¬ 
volves eight vibrational hg modes, i.e., 40 vibrational 
coordinates.^ The corresponding vibronic coupling 
parameters for have been recently extracted from 
DFT calculation,l4Sl while the reliability of this approach 
was proven by a satisfactory reproduction of photoemis¬ 
sion spectrum for C^q. 14^ Nevertheless, in the present 
calculations the use of a full multimode description of 
JTE on fullerene sites seems to be impractical. For 
this reason, the eight-mode JT interaction on fullerenes 
has been replaced with an effective single-mode one 
©■ Thus the two parameters, hw = 87.7 meV and 
g = 1.07 were obtained via the reproduction of the 
JT stabilization energy and the energies of the lowest 
vibronic excitation of ion.^^l In the model JT 

Hamiltonian the quadratic vibronic couplings are 
not included because, as we discussed in Ref. HOl they 
are weak in Ceo anions and do not give significant effect 
on the JT dynamics of in cubic fullerides. Finally, 
the Hand’s rule coupling parameter, Jr = 44 meV, was 
also taken from the DFT calculations .1451 This is not the 
case of interelectron repulsion parameters of fullerene 
site, which are strongly renormalized by screening in 
fullerides .mi In the present work, the Coulomb repulsion 
U is treated as a free parameter. U is defined here as 
the average repulsion of two electrons in Cg^ for a cubic 
(undistorted) LUMO band: 


[/= I ([/|l+4[/i) + (5) 
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A. Adiabatic orbitals 


The hg normal coordinates on each site m, ^^ 7 , are 
expressed by polar coordinates, (gm, 7 m, ^m, 0 m )-^ 
Introducing a unitary matrix, 

Sixi^m) = [Bp{-fr^)Cp{Or^)Dp{(t)^)],^ , ( 6 ) 

we transform the electronic basis (A = x, y, z) into adia¬ 
batic basis (I = 1, 2, 3) on each Ceo site (Fig. |^, 

~ ^ ^^A(^m)c|^Acr* C^) 

X=x,y,z 


Here, = (7m, ^m, 0m), and Bp^ Cp, Dp are the Euler 
rotational matrices defined in Ref. [271 By the transfor¬ 
mation of the electronic basis, Eq. Q, the linear vibronic 
term of the JT Hamiltonian 0 becomes diagonal: 


Uljt = S^UljtS 

= EE Hujgqn 


cos (am + 


^mlcr 


+ COS 


(“” 3 ) 


I rini2cr COS Q^m^mScr 


(8) 


where S = Hm and Sm is the unitary operator whose 
matrix element is given by Eq. (§. Eq. ^ shows that 
the amplitude of the JT distortion is determined by ra¬ 
dial coordinates and o^m, and the direction of the JT 
distortion in the space of the five dimensional hg normal 
coordinates is defined by Euler angular coordinates 0^ 
(Fig. §. In the described coordinate system, the elastic 
energy term in Eq. ^ is written as 

Uel = E 

m 

i.e., is invariant under the unitary transformation (§. On 
the other hand, the kinetic energy term changes, which 
is discussed in Sec. HYU 

Under the transformation of the electronic basis 0 , 
the transfer Hamiltonian ([^ becomes 

k = S^HtS (10) 

m,m' ll'a 

where is 

XX'=x,y,z 

( 11 ) 


We note also that ^bi is invariant under the unitary 
transformation ^ due to the isomorphism of LUMO 
shell of Cgo" to the atomic shell. 

Eor any Euler angles, (2^, the JT potential term, Eqs. 
Q and (|^, and the bielectronic term has the same form. 
Therefore, the adiabatic potential energy sur face o f an 
isolated Cgg" has continuous minima (trough)^^!!] 
in the presence of the term splitting. In the case of , 



FIG. 3. (color online) Ordered fee KsCeo- The orange ball is 
Ceo, and the red and blue spheres are K atoms in octahedral 
and tetrahedral interstices. 


the potential surface has three dimensional (3D) trough 
at 

and a = 1^12. Substituting cj, and Jh above into 
Eq. ( p!^ , the amplitude of the JT distortion is g = 
0.989 X \/3g, indicating that the effect of the Hand’s rule 
coupling on the JT potential surface of Cgg is small. 


III. GUTZWILLER APPROACH TO STATIC 
JAHN-TELLER SYSTEMS 


A. Self-consistent Gutzwiller approach for the 
LUMO bands 


The merohedral disorder in the K 3 C 60 lattice and the 
orientation of the JT distortions on the fnllere ne site s do 
not have important effect on the band energy.1^^1^ The 
change in Hartree-Eock energy per Ceo site due to the 
disorders is only 14 meV,l^ which is smaller than the JT 
energy of CgQ by one order of magnitude. The varia¬ 
tion will be further reduced by the electron correlation 
as is discussed in Sec. |!vb| Therefore, for the sake of 
simplicity, we further consider a K 3 C 60 in an ordered fee 
lattice (Eig. [^. As a possible scenario of static JT effect 
we consider equal JT distortions on fnllerene sites of the 
following form: 

(gm, 7m5 0rn) ~ (^, 0, 0, 0) , 

or in conventional coordinates: 

(^m0, ^me, ^m77, ^m(’) — (0, V^g, 0, 0, 0), (U) 

which do not remove translational symmetry of the lat¬ 
tice. Here, g is a variable. The direction of the JT dis¬ 


tortion (13) corresponds to the one which gives the max¬ 


imal static JT stabilization in the case of isolated C; 


i3- 

'60 
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ion JHEH Under the distortion (13), the adiabatic orbitals 
I = 1, 2,3 correspond to x, z, respectively (Fig. [^, and 
the linear JT term dsl) reduces to 


b^LJT = EE («!] 


^mya 


)• (15) 


This expression shows that the 2 : orbital level remains 
unchanged while the x a nd the y levels are stabilized and 
destabilized, respectively! ^^ * ^^ * 

The Gutzwiller wave function, |Tg), is expressed as 


|^g) = Pg|^s), 


(16) 


where |4>s) is a Slater determinant, and Pq is a 
Gutzwiller projector. The Slater determinant is written 
as follows: 


i^s)=n4<.io)’ 

pkcr 

oikm 


=E 

mA 


'^pkcr 


'^A,pkc|^Acr’ 


^/N 


(17) 

(18) 


where p is a band index, N is the number of sites in the 
system, and ux^pi^ is a variational orbital coefficient. We 
note that the band described by |4>s) is not constrained to 
obey the cubic symmetry. In order to include properly 
the effect of JT distortions and of electron correlation, 
the variational parameters (^4) in Pq have to be orbital- 
specific: 


-Pg = P exp 



^AA'f^mAcr'hmA'cr' I 5 (^9) 


where ^aa' are real and symmetric with respect to 
interchange of indices. Therefore, the projector (19) 


is described by six independent Gutzwiller parameters 
{Axx') instead of a single param eter used in conventional 
Gutzwiller wave fnnction. h^FQ I a general case. A, A' 
denote natural orbitals on the site m. For the chosen 


JT distortions (13), preserving the orthorhombic site- 
symmetry, these natural orbitals coincide with the or¬ 
thorhombic X, p, 2 : tiu LUMO orbitals. Due to equal dis¬ 
tortions (13) on all fullerene sites, Axx' are independent 


from the index m. 

The calculations of expectation values with |Tg) have 
been done within the Gutzwiller’s approximation.^^^^^ 
Within this approximation, the energy per site. 




1 (TgI^I^g) 

N (TgI^g) ’ 
consists of the band energy, 

Et = ^ Qxx'Txxq 
XX'a 

the elastic energy the linear vibronic energy, 

Gjt = '^-^fuvgq (ux - riy), 


( 20 ) 

( 21 ) 

( 22 ) 


and the bielectronic energy E\^[. Here, qxx' is the 
Gutzwiller’s reduction factor, rxx' is 

1 


^aa' — ^^^aa'Paag 


where is the Fourier transform of 


- E 


—ik-Am^Am 
^ ^AA' 5 


Am 


is the density matrix at a k point. 


occ 

PXX' = ^'^A,pk'^A',pk, 
P 


and nx (A = x, p, z) is the occupation number, 
(^G|fimAcr|^G) 


nx = 


(23) 


(24) 


(25) 


(26) 


(^gI^g) 

The explicit forms of the occupation number nx, the 
Gutzwiller’s reduction factor qxx', and the bielectronic 
energy F’bi are given in Appendix The Gutzwiller 
projector does not influence the on-site density matrix, 
hence, Eq. (I^ corresponds to 


(27) 


n\ = (^slnmAal^s) = 


Hereafter, we use the form (27) for nx- 


The ground state for different amplitudes of JT distor¬ 
tion q is obtained by minimizing the energy per site (20) 


with respect to {ux,p\^} and {Axx'}, which is performed 
in two steps. The first one is the variational calculation 
of Eg with respect to {ux,p\i} for fixed {Axx'}- The re¬ 
sulting self-consistent equations in the case of static JT 
effect are obtained in the form: 




AA'^A',pk — ^pk'^Ajpk: 


(28) 


where the one-particle Hamiltonian is 

dq, 


hx\' = qxx'txx' + hx' 


E 


dn 


KK' 1 dE]ji 

+ t:' 


2 dn\ 


- ^Hujgq{6xx - hy) 


(29) 


and e^k is the Gutzwiller’s orbital energy. Using the so¬ 
lutions of Eq. ( [ 2 ^ , {ux^pi^}, the occupation numbers 
{nA} are recalculated via Eq. ( [27| ). The chemical poten¬ 
tial is found by consecutive population of Gutzwiller’s 
orbitals following the aufbau principle. The second step 
is the minimization of Eg with respect to { Aaa' } for fixed 
{^A,pk} and {nx}, 


dE, 

dAxx 


^-E 

(JKK‘ 


dq^ 


dAxx' 


dEu 

dAxx' 


= 0 , 


(30) 


using the numerical algorithm proposed in Ref. |311 The 
two minimizations, (28) and (30), are repeated itera¬ 


tively until variations in the occupation numbers and the 
ground state energy become smaller than thresholds. 
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FIG. 4. (color online) (a) Total energy Eg{q) of KsCeo as 
function of amplitude of static JT distortion (13) for several 
values of U § . The dashed line indicates the minimum of 
Eg. (b) Energy components contributing to the total energy. 
UhjT is the linear JT energy ^ and Ue\ is the elastic energy 
^ in ^jt; Et is the band energy (21), and E^i is the bi- 
electronic energy, is the ground state adiabatic potential 
of an isolated Cg^. Eg and E^i are set to zero at Q— O- T he 
unit of q is the amplitude of zero-vibration of hg mode.^^ (c) 
Occupation numbers of LUMO orbitals nx as function of q. 
X, y, and z are orbital components. 


B. Static Jahn-Teller instability in KaCe 


The ground state energy Eg ( [2Q|) a s a function of the 
JT distortion q is plotted in Fig. ]^. Quite unexpect¬ 
edly, the energy curve Eg{q) has two minima, one at the 
undistorted configuration q = 0 and the other at g ~ x/3g 
(= 1.85) which corresponds approximately to the am¬ 
plitude of JT distortion in an isolated Cgg ioii- For U 
smaller than the critical value Uc = 670 meV the static 
JT distortion is quenched (g = 0). The minimum cor¬ 
responding to JT-distorted sites lowers with the increase 
of [/, at = Uc the values of the two minima equalize 
and for U > Uc the JT distortion achieves its equilib¬ 
rium value matching approximately the distortion in an 
isolated Cgg • 

The static JT effect has been investigated in A^Cqq 
within the local density approximation (LDA) of DFT, 
for which completely quenched JT distortions have been 
found.l^ Since the static JT effect in C 


gQ is stronger than 


in CgQ anion,!^ it was concluded that the JT distortions 
in AsCeo are also quenched. However it was recently re¬ 
vealed that the LDA calculations underestimate the JT 
stabilization energy of C^q by ca 30 On the other 
hand the broken-symmetry Hartree-Fock (HF) calcula¬ 
tions predict smaller Uc for static JT instability than the 
present calculations and orbital dispr oport ionation of the 
intrasite charge density in fullerides reason of 
this discrepancy is that the broken-symmetry HF calcula¬ 
tions exaggerate the tendency towards the stabilization 
of low-symmetry electronic phases. Thus the splitting 
of the LUMO band is overestimated and is mainly con¬ 
tributed by the interelectron repulsion,!^ suggesting that 
an approach based on a single Slater determinant is not 
flexible enough to include properly the effects of elec tron 
correlation in orbitally degenerate bands (see Sec. V A 
for detailed discussion). 

The structure of Eg{q) is mainly determined by the 
g-dependent contributions, the band energy Et (21) and 


the JT potential Uljt which are found in competi¬ 
tion (Fig. 1^). The contribution of the band energy term 
is the largest when the orbitals are hybridized and equally 
populated, which takes place in the weak correlation limit 
{U ^0). On the other hand, the contribution from the 
JT term is the largest when the disproportionation of the 
electronic charge among the LUMO subbands, accompa¬ 
nying the JT distortions, is full, {rix^ny^nz) = (1,0,1/2). 
At the same time, the splitting of the orbital levels pre¬ 
vents the hybridization and vice versa. As Fig. & sug¬ 
gests, the ground state energy Eg{q) consists of two po¬ 
tential energy surfaces which cross at g = gc- For small 
distortions q < Qc^ the band energy exceeds the JT en¬ 
ergy and the latter is quenched compared to the case of 
an isolated (^id)’ while for g > gc the JT energy 
takes over. Because of the hybridization, the occupation 
numbers nx are fractional (Fig. [^) and the linear JT 
energy Uljt is not proportional to g unlike the isolated 


CgQ molecule. The band energy is reduced by the intra¬ 
site Coulomb repulsion by quenching charge fluctuations 
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on Ceo’s. Since the band energy in RbsCeo and CssCeo 
is smaller than in KsCeo, while U is larger pi the static 
JT instability is favored even more in these fullerides. 


IV. GUTZWILLER APPROACH TO 
DYNAMICAL JAHN-TELLER SYSTEMS 


A. Dynamical Jahn-Teller contribution 


Another important ingredient is the energy gain aris¬ 
ing from the dynamical delocalization of JT distortions 


at each Cgg anion. The energy gain in isolated C; 


3- 

60 


amounts to ca 90 meV which is more than a half of 
the static JT stabilization (ca 150 meV) in this anion.^ 
To assess this energy gain in fullerides, one should take 
into account that the JT effect on Cgg sites in fullerides 
is different from the case of isolated fullerene anions. 
The main difference is that the LUMO orbitals on the 
fullerene sites do not have the same populations as in 
an isolated C^q (Fig. [^), which leads, in particular, 
to lower values of the amplitude of dynamical JT de¬ 
formation in fullerides. Only in the case of full orbital 
disproportionation (Sec. IIIB) the deformation achieves 
the equilibrium value in a free ion {q = 1.85) and the cor¬ 
responding energy gain owing to dynamical JT effect is 
maximal. One should stress that in the case of dynamic 
JT effect the adiabatic orbitals I = 1,2,3 on fullerene 
sites are not fixed electronic orbitals A = consid¬ 

ered in the previous section but are their linear combi¬ 
nations with Ojn-dependent coefficients, Eq. 0 Pto 
simulate the dependence of dynamical JT effect on the 
extent of orbital disproportionation, we introduce the ef¬ 
fective vibronic coupling constant. 


geS = gini-n2), (31) 


which varies from 0 to ^ when the orbital disproportion¬ 
ation Til — 77-2 varies from the minimal value (0) to the 
maximal value (1). Note that ns has an unchanged value 

1 / 2 . 

Diagonalizing the JT Hamiltonian ^ for different val¬ 
ues of and extracting the ground state energy at 
corresponding static JT distortion, ^fiujgl^/2^ together 
with the energy of zero-vibrations at distorted point, 
bflujl2^ we obtain the dynamical contribution, Edjt, to 
the ground vibronic level. Figurej^shows the dependence 
of this contribution on ^eff for the case of effective single¬ 
mode JT Hamiltonian of Cg^.^Note that the existence 
of the energy gain due to dynamical delocalization of JT 
deformations does not guarantee by itself the develop¬ 
ment of dynamical JT effect on fullerene sites. For the 
latter to take place, an additional condition are the small 
variations of the band energy under arbitrary JT distor¬ 
tion on Cg^ sites, which is investigated below. 



FIG. 5. (color online) The gain of JT stabilization energy 
(meV) due to dynamical JTF in function of ^eff- 


B. Form of the ground vibronic state 


Previous ab initio investigations have shown that 
the low-lying vibronic states in an isolated Cgg 
can be described satisfactorily within the adiabatic 
approximation.!^ This approximation can be extended 
over the AsCeo crystal. Following the molecular 
Mproach,!^ first we perform the unitary transfor mation 
Q to diagonalize the linear vibronic term in Eq. 


k = S^HS = 4ad + 4ot + + i?bi, 


fj _ ^ 

^rad — 2_^ ^ 


^-4 


+ 


' dqrr 
d 


^el 

d 


(32) 


sin 


sin 3ar, 


d 

dan 


Huj 9 


Hrot — 


+ 








+ 




^ sin^((am — 27r/3) sin^((am H-27r/3) 


sin^ 


(33) 


Here, is the sum of the linear vibronic term (|8 


and the transfer part (10), nuclear 

angular momenta in the initial orbital basis (A^,, A^, in 


Ref. |27l respectively), and L^j{j = 1, 2, 3) are electronic 
angular momenta: 


fel _ 


rel _ 
^m2 — 


fel _ 
^m3 — 


^ (^i!n2cr^m3cr c|^3cr^m2cr^ 5 

(j 

(34) 

^ (^^m3cr^mlcr “ 5 

(7 

(35) 

^ (^mlcr^m2cr “ ^ln2a^rnla^ • 

(36) 


(T 


For arbitrary JT deformations on sites, the system does 
not possess translational symmetry anymore. In the case 
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of intermediate to strong vibronic coupling, the ampli¬ 
tude of dynamical JT deformation go is not small. Since 
the LUMO orbitals of each fullerene are, on average, oc¬ 
cupied by three electrons, the vibronic term has a mini¬ 
mum at Of = 7r/2.l^ Substituting 


qin = qo + q'm, 


dm — — + 


(37) 


into Eq. (32), we obtain 


^ = 4ad+i7rot+^e7+^l 
d 


TT fUjJ 

-Hrad = 2^-^ 


bi? 




dq'm V ”^ 9 ; 


d 

m 


+ 


1 52 

ql da'^ 


o V—' /2 

N^qi + -yC, 


= +4(r™i+r«j,2) 


(38) 

qm=qo 

(39) 
2 


47 = EE- ^hbjgqo {hmla - nm2a) 

m a 

+ E 

m,m' ll'(T 


(40) 


(41) 


where g^ and are the deviations from the equilibrium 
point. The above derivation is based on the assumption 
that the radial JT coordinates (37) remain unchanged 
under the electron transfer. The justification for that, 
for the n eglect of JT polaronic effect will be given 
Following the adiabatic approximation,!!! 


i.e., 


in Sec. VIC 


m 


Eq. (38) the terms smaller than l/gg are neglected and, 
consequently, the radial degrees of freedom (g', a') are de¬ 
coupled from the other degrees of freedoms correspo nding 
to the rot ation of JT deformation in the 3D trougli^!!^ 
(see Sec. IIA for the trough). The rotational Hamilto¬ 
nian (|40|) has nonadiabatic terms, V = fiujL^^L^-/qQ. 


Neglecting these terms,!!! we obtain 


Hp,r\ - Hy- 


r(i) 




nuc \ 


I rVnuc I rVel i tt 

T T HyQi T 

m L 


(42) 

(43) 

(44) 


where + L ^2 + The adiabatic approxi¬ 

mation is valid when the energy gap between the ground 
and the first excited energies AE is large compared with 
the matrix element of the nonadiabatic term \V\. The 
ratio of \V\ ^ fujul{3g‘^) and AE ^ 3flujg‘^/2 for is 
|V|/AE ^ 1/5, which justifies the application of adia¬ 
batic approximation in the present case. In this estima¬ 
tion, a value \L\ ^ 1 was taken. 


Diagonalizing (|41|), the Hamiltonian is written in 
the basis of adiabatic 


Dand orbitals: 


h: 


(1) 




(45) 


where D = {Dm} is the set of all Euler angles on all Ceo 
sites in the system (Fig. [^, i indicates adiabatic band 
orbital, denotes its energy, and is given by 




mla ’ 


(46) 


I 


For the ordered system (13), the coefficient Umii reduces 
to e'^^'^ux^p\^/^/N appearing in Eq. (18). 


Then the solution of the Hamiltonian (42) in the adi¬ 


abatic approximation for the ground and low-lying vi¬ 
bronic states has the form: 

\^iR, n)) = (47) 

where |4>g^) is the Slater determinant of occupied adia¬ 
batic band orbitals (46): 


i^s 7 ^^))=n 4 (^^)io), 


(48) 


and and x^°^(D) are nuclear wave functions de¬ 

pending on radial R = {g(^, a'^} and rotational D nuclear 
coordinates, respectively. The factorization of nuclear 
wave function became possible due to the separation of 
radial and rotational degrees of freedom in the adiabatic 
Hamiltonian (42). Furthermore, the radial coordinates of 


different sites are independent from each other (see Eq. 
(39)), hence, the radial part is the product of the 
ground vibrational wave functions of all sites: 




(49) 


Further calculations are greatly simplified under the 
assumption that the dependence of (46) on Euler 


angles is relatively weak. This seems to be the case when 
correlation effects become important, leading to signifi¬ 
cant reduction of band energy (Fig. |^) and strong 
separation of Gutzwiller bands (x, ^ in Fig. for ho¬ 

mogeneous JT distortions (p!^). Indeed, the hybridiza¬ 
tion of the adiabatic orbitals in this case mainly arises via 
resonant interactions (Fig. because the width of in¬ 
dividual bands is small compared with their Jahn-Teller 
splitting (150 meV), thus resulting in a weak mixing of 
the off-resonant adiabatic orbitals. The hybridization 
arising from resonant interactions should be weakly de¬ 
pendent on transfer parameters. For example, in the case 
of two-site model (Fig. [^, the “band” splitting of pairs 
of interacting resonant adiabatic orbitals Ae is strongly 
dependent on the Euler angles on two sites, 

A€^\24^{nA,^B)\, (50) 

while the adiabatic “band” orbitals. 


|'0icr(D)) ~ ^ T 1^) 


(51) 
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FIG. 6. Hybridization of adiabatic orbitals in a two-site 
model. 


have angle-independent mixing coefficients. 

Neglecting the (^-dependence of coefficients Ur^u in Eq. 
(46), the eigenvalue problem for the pseudorotational nu¬ 
clear wave function reduces to the equation: 


' T^nuc I 77'el I ttx 




where Eq^{Q) is the adiabatic band energy: 


(53) 

i(j 


and is the expectation value of ^4* (44): 




kad\ 


(54) 


The direct calculation of this matrix element gives: 

^2 

-^rot ~ ^ ^ ^ 2 10^m2 IhUjnS 

m ^^0 

- 4nminm2 - 16nm2^^m3 “ 16nm3^^ml) 

occ ^2 

+ EEi(4|t^mu|^|C^2i|^ 

m i 

+ 16 |C/^2i|' |C^3i|" + 16 |C^3i|' ,(55) 

where Um/ are populations of the adiabatic orbitals (/) on 
the site m. The last term is smaller than the other terms 
by l/iV because \U^u\ ~ l/\/iVandEi « 

l/A/", while the occupation number Um/ = ~ 

1. Neglecting the last term, we obtain 


of electron density among three orbitals, (ni, 712 ,^ 3 ) = 
( 1 , 0 , 1 / 2 ) (Sec. [m^ Note the lack of (l-dependence of 
the energy in Eq. (56), which is the result of neglected 


0-dependence of the coefficients Umii ia Eq. (46) 


On the other hand, the adiabatic band energy Ef (53) 
is 0-dependent even if the coefficients Umii are not, and 
this dependence a priori is not weak. This 0-dependence 
is estimated here by direct calculations of the uncorre¬ 
lated band energy E^ for different directions of ordered 
JT distortions, q = V^g. The obtained variations of E^ 
do not exceed 12 meV (Eig. [^. The variation of E^ 
will be even smaller for disordered system because the 
Euler angle dependence is smeared out by the disorder. 
Including electron correlatio n effec ts via the Gutzwiller’s 
ansatz described above (Sec. |mB| will result in the case 
oi U = Uc (corresponding to q = 1.85, see Eig. Ef) to 
a reduction of uncorrelated E^ —240 meV) by one 
order of magnitude (Eig. Ef). At the same extent will 
reduce the variations of the band energy in function of 
the direction of JT distortions, which means that they 
are negligible compared to the dynamical contribution 
to JT stabilization energy (Eig. |^. 

In A 3 C 60 crystals the JT pseudorotations after the Eu¬ 
ler angles Q can be also hindered by intermolecular vibra¬ 
tions. However, the energy of these vibrations (?^ 5 — 10 
me VES) 

is much lower than the energy gain due to delo¬ 
calization of JT deformations in the trough. 

Hence the vibronic dynamics is expected to be un- 
quenched, like in insulating fullerides Cs 3 C 6 of^ Given 
the near independence of the band energy Eq^ on the 
pseudorotation coordinates of sites ((1), and the full 
(1-independence of the contribution ( [5^ , the pseudorota¬ 
tional Hamiltonian (52) becomes merely a sum of on-site 


contributions. Each such contribution is an operator de¬ 
pending on dm Euler coordinates of the corresponding 
site, Eq. (43), which means that the pseudorotational 


wave function factorizes. 




(57) 


with being eigenfunctions of one-site operators in 
Eq. ( [ 4 ^ . Then, taking into account the factorization of 
the radial part, Eq. (49), the Gutzwiller wave function 


with dynamical JT effect on fullerene sites has the form: 


|v|/G)=PG|^r)xnXm 


(58) 


E 


el 

rot 


E 


2 “1“ 107T/ni2 “1“ 

8^0 


471rnl7^m2 • 


(56) 


and the Gutzwiller projector (19) will involve now pop¬ 
ulation operators for adiabatic orbitals on the fullerene 
sites: 


The obtained energy is additive over the sites, with one- 
site contributions being equivalent with the correspond¬ 
ing energy of an isolated Cg^, = 5^^/(4gQ) (Eq. 

(37) in Ref. [27|) . in the case of full disproportionation 


-Pg = H 

m 



E 


A. 


dl' a' 


(59) 
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FIG. 7. (color online) Dependence of uncorrelated band 
energy on the Euler angles of JT distortions 0 (a) and 7 
(b), respectively, for q = \/3g = 1.85. 


and ^djt ^be dynamical JT contribution 


TT'ad _ 

^DJT — ~ 


3Huj 


3fujj 

W 


■ Er-^ 


(62) 


eff 


The zero-point energy of the five-dimensional harmonic 
oscillator is set to zero. The first and the second terms in 
Eq. appear from the radial Hamiltonian (39 and 


Ej-ot is the eigenvalue of the pseudorotational Hamilto¬ 
nian ( [ 5 ^ . Furthermore, the dynamical contribution (62) 
is replaced by the exact F^djt (Fig. |^, yielding 

E = Et + Eu- + ^DJT(5eff). (63) 

The ground state for dynamical JT system is obtained 
by self-consistent minimization of the energy (63) with 


respect toJux,p\i} and {Hyy/}. We obtain similar formula 
as Eqs. (28) and (30), with the only difference in the JT 


term of one-particle Hamiltonian: 


hxx' = qxx'txx'^^xx 


dq^ 


IdE^ 


bi 


3fwj 

+ I-+ 2 


^ dnx '2 dnx 

. kk' 

g 5F^DJT(^eff)^ 


%eff ) 


(64) 


D. Dynamical Jahn-Teller instability in KaCeo 


C. Self-consistent Gutzwiller approach for the 
ground vibronic state 


The ground state energy of the dynamical JT system 
is obtained by minimizing the total energy per site. Al¬ 
though the adiabatic band orbitals (46) correspond to a 


disordered system, this will not pose any complication if 
we assume that the band energy of these orbitals is in¬ 
dependent on the form of adiabatic orbitals, i.e., on the 
three Euler angles characterizing the “direction” of JT 
distortions on sites. This seems to be indeed the case 
given the weak dependence of band energy on the local 
JT distortions established above (Fig. [^. Then the cal¬ 
culation of the electronic part of the energy can be done 
for a particular case of Euler angles equal on all sites, 
yielding the previous result for a translational system, 
while the nuclear part of the wave function (58) will give 
the dynamical contribution. 


Hence, within the adiabatic approximation (42), the 


ground energy with the Gutzwiller’s wave function (58) 
is given by 


F/ad — F/t T F/bi 


Sfkog^g 


■E\ 


ad 

DJT5 


(60) 


Minimizing the total energy (63), we obtain the ground 


energy in the presence of the JT dynamics on sites (Fig. 
1^). In the case of static JT effect, the JT distortion ap¬ 
pears for U > 670 meV (Fig. |^). We can see, however, 
that the JT dynamics enhances the dynamical JT defor¬ 
mation, and, consequently, the disproportionation of the 
occupation numbers in the adiabatic orbitals are also en¬ 
hanced (Fig. 1^). As a result the critical value of electron 
repulsion parameter for JT instability {Uc) is significantly 
reduced in the dynamical case. In particular, the criti¬ 
cal value is smaller than the estimated 7/ = 750 meV for 
KsCeoP^ hence, the metallic fullerides always exhibit dy¬ 
namical JT instability in the ground state. This explains 
the absence of staggered JT deformations in the x-ray 
diffraction data of AsCeo- Furthermore, since U > Uc 
the equilibrium JT distortions on sites will be close to 
maximal possible, i.e., to their values in a free Cgg 


^ ion. 


V. EFFECT OF ELECTRON CORRELATION 
AND JAHN-TELLER INSTABILITY ON 
ONE-PARTICLE STATES 

A. Orbital disproportionation 


where the dynamical JT deformation qq is replaced by 
go = VSges, (61) 


The electron correlation and the JT effect induce dif¬ 
ferences in the population of the three LUMO orbitals 
on fullerene sites (orbital disproportionation). Within 
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(a) 




FIG. 8. (color online) (a) Total energy and (b) occupation 
numbers of LUMO orbitals nx with static (blue) and dynam¬ 
ical (red) JT effects as a function of U. The blue points and 
the dashed line correspond to the minimum with static JT 
distortion and the global minimum of Eg, respectively. Total 
energy at g = 0 is set to zero at each U. x, y, and z in (b) 
are orbital components under distortion (13). 


the broken-symmetry Hartree-Fock approachp^ the JT 
and the bielectronic energy per site is 








(65) 


where n is the total population of the fullerene site, An/ 
is the deviation of the occupation of the orbital subband 
from the case of cubic symmetry (1/2), and one single 
average electron repulsion parameter U (§ is used for 
simplicity. The HF energy with full disproportionation, 
(ni,n2,n3) = (1,0,1/2), is lower than the energy of the 
degenerate system (ni,n2,n3) = (1/2,1/2,1/2) by 


= - 


3Hujg‘^ 


U 


( 66 ) 


The orbital disproportionation is seen also in the present 
Gutzwiller’s treatment (Fig. [^. In terms of the elec¬ 
tron configurations, the equal population of three LUMO 


bands in a cubic band structure results in their equal 
probability (1/2^). The HF type symmetry breaking 
equally enhances the weights of four configurations, ipf, 
(both spin projections) and and quenches the 

others, leading to the gain of bielectronic energy per site 
^bi by U/2. The weights of configurations among 

the four are further enhanced and the rest of them are 
further reduced in the Gutzwiller treatment, which ad¬ 
ditionally lowers E]ji by 17/4 in the limit of strong cor¬ 
relation. The latter becomes possible because of multi 
determinantal character of the Gutzwiller ansatz. 

Despite the larger gain of E\y[ in Gutzwiller approach 
compared to HF one, the latter predicts smaller Uc for 
the static JT distortion. This is due to the artifactual 
feature of the broken-symmetry HF approach mentioned 
above which leads, in particular, to orbital disproportion 
in K3C60 without JT effect on fullerene sites.l^ Indeed, 
even in the absence of the vibronic coupling, ^ = 0, the 
broken-symmetry HF state is more stable than the cubic 
band solution by U/2, Eq. (66). On the other hand, the 
Gutzwiller’s wave function is not disproportionated in the 
absence of JT effect, which is testified by equal popula¬ 
tion of three LUMO orbitals at q = 0 point for arbitrary 
U (Fig. 15). This is the result of a higher flexibility of 
the Gutzwiller’s wave function, which can include various 
configurations without changing the bielectronic energy, 
such as equally populated configurations of type. 

Orbital disproportionation can be directly observed in 
spectroscopy, e.g., in photoemission spectra of fullerides. 
Following the preceding discussion, the quasiparticles will 
belong to subbands with definite orbital index, I = 1, 2, 3, 
separated by energy gaps (Fig. 15). The centers of grav¬ 
ity of these subbands is expected to coincide with the 
centers of Gutzwiller subbands obtained as solutions of 


Eq. (28). The latter are expressed by the sum of the JT 


splitting and the Coulomb repulsion energy: 


ebi,z = ^ U 




(67) 


where is the ground state wave function. Consequently, 
the energy gap between centers of weight of the subbands 
is expressed as: 


^^disp — 




eff 


Aebi. 


( 68 ) 


The bielectronic part Aebi of Eq. (68) for br oken- 
symmetry HE solution is given by Ae^ = U/2^^^ 
Aebi for Gutzwiller wave function is calculated using Eq. 
(B5). Ae^^ and Aebi for Gutzwiller’s wave function are 
shown in Eig. Ae^^ monotonically increases with 

U, while Aebi for the Gutzwiller’s solution approaches 
to zero. The bielectronic contribution Aebi becomes zero 
because the system approaches to the isolated molecular 
limit: when electrons are completely localized due to the 
metal-insulator transition, the splitting of the subbands 
reduces to the JT splitting in isolated Cgg We can 

see from Eig. 15 that Aebi, while exaggerated in broken- 
symmetry HE approach, is not an artifactual feature but. 
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(a) (b) 



FIG. 9. (color online) (a) The energy gap between Gutzwiller 
subbands, Acdisp (68). (b) The bielectronic part Acbi of Eq. 
( [6^ as function of U (meV) for HE (dashed line), Gutzwiller 
with dynamical JT (solid line), and Gutzwiller solution with 
static JT distortion (points). Acbi for the global minimum 
with the static JT effect is shown by the red dashed line. 



on the contrary, gives a non-negligible contribution to the 
splitting of quasiparticle subbands in the metallic phase. 
Figure 10 shows that the charge fluctuations (probabili¬ 
ties of configurations with n = 2,4) is suppressed at [/ < 
700 meV, signaling the arising of metal-insulator transi¬ 
tion. 


B. Density of states of uncorrelated LUMO band 

It is also of interest to find out how the uncorrelated 
band structure is affected by JT instability. Figure [TT| 
shows the density of states (DOS) of the uncorrelated 
LUMO band in the presence of equilibrium JT distortion 
{q = 1.85). Compared to cubic band structure, we see 
a strong enlargement of the bandwidth by ca 300 meV. 
The analysis of partial density of states shows that the 
degenerate LUMO band splits into three subbands (Fig. 
pJ^) mainly contributed by one of the adiabatic orbitals 
(these are x, y and 2: for the distortion @). This means 
that the electron correlation in fullerides does not take 
place in a degenerate LUMO band. In particular, the 
Mott-Hubbard transition in cubic fullerides basically oc¬ 
curs in a split band structure, where half-filled is only 
the middle band. This calls for reconsideration of the 
role played by orbital degeneracy in the Mott-Hubbard 
transition in fullerides. 


VI. DISCUSSION AND CONCLUSIONS 

The vibronic interaction and the electron correlation in 
^sCeo are concomitantly treated by a new approach pro¬ 
posed here based on self-consistent Gutzwiller’s ansatz 
with orbital-specific variational parameters. The present 


FIG. 10. (color online) The p robab ilities of the electron con¬ 
figurations appearing in Tq, ^y p^ l ^^ las functions of U. Dotted, 
dashed and dot-dashed lines correspond to cubic symmetry 
band structure, solid lines correspond to dynamic JT effect, 
and symbols correspond to static JT effect. Black, blue, and 
red indicate the electron configuration ' 0 i'p 2 ? sum of z/’s over 2 
(or 4) electron configurations, and sum of z/’s over 3 electron 
configurations, respectively. The vertical dashed line indi¬ 
cates Uc for the static JT instability, z/’s for 2 and 4 electrons 
are almost identical to each other (only the data for n — 2 
are shown), z/’s for n = 0,1, 5, 6 are not shown here because 
they are close to zero. 


Gutzwiller’s calculations with realistic vibronic con¬ 
stants, Hand’s rule coupling and parameters of the 
LUMO band predict that both the static and the dy¬ 
namical JT deformations arise in AsCeo- Since the elec¬ 
tron correlation quenches the band energy, the localiza¬ 
tion of the electrons is enhanced, and consequently, the 
JT distortions on Cgg sites is facilitated. It is shown that 
the dynamical JT instability appears for smaller on-site 
Coulomb repulsion, U < 500 meV than the static one 
(Fig. [^). Due to the existence of the dynamical JT 
distortion, the adiabatic LUMO band splits into three 
subbands (Fig. pff] ). An indirect experimental evidence 
for the existence of dynamical JT effect in fullerides is 
given by NMR spectroscopy of CssCeo, showing that fea¬ 
tures attributed to dynamical Jahn-Teller effect in its in¬ 
sulating phase persist when this material is brought into 
metallic phase by applying an external pressure.^ 
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(a) 




FIG. 11. (color online) (a) DOS’s per Ceo and spin for the un¬ 
correlated LUMO band of KaCeo in the absence (dashed) and 
the presence (solid) of equilibrium JT distortions on fullerene 
sites, (b) DOS with the JT distortion and partial DOS’s 
corresponding to three adiabatic orbitals. The vertical lines 
indicate Fermi levels. 


(a) 



U (meV) 


(b) 



A. Correlation in split bands 

The results of the present work do not support the es¬ 
tablished view that the electron correlation in fullerides 
takes place in a degenerate LUMO band. As was shown 
by Gunnarsson et a/.j^LrP^and Han et al^ the orbital 
degeneracy of the band/metal sites increases the critical 
ratio U/w for Mott-Hubbard metal-insulator transition 
where w is the width of the band. Gunnarsson et al. 
has found that this ratio is 1.5 — 2.5 for A3 Ceo, which 
is significantly larger than the critical ratio U/w ^ 1 for 
Mott-Hubbard transi tion in lattices with orbitally non¬ 
degenerate sitesWith the bandwidth w ~ 0.5 — 0.6 
evEzl (Fig. [TT^ ) and the estimated U ^ 1.3 — 1.6 e\®l 
it was natural to conclude that the orbital degeneracy of 
the LUMO band, leading to large critical values of U/w, 
is the reason for KsCeo and RbsCeo to remain metals.^ 
This picture has become a basis for the interpretation of 


FIG. 12. (a) The probabilities of the electron conhgura- 

tions appeari ng in Gutzwiller wave function of half-hlled cu¬ 
bic system, ^yliniMl functions of U. The dashed lines and 
solid lines correspond to Us for non-degenerate (one electron 
per site) and three-fold degenerate systems (three electrons 
per site), respectively, (b) The DOS’s for the non-degenerate 
(solid line) and degenerate (dashed line) systems. Vertical 
lines show the position of the corresponding Fermi levels. 


metal-insulator transition in fulleridesj ^^*^^ * in particular, 
in Cs3 CeoContrary to that, the JT-split correlated 
state derived here exhibits the Mott-Hubbard transition 
at a lower critical ratio U/w. Indeed, Fig. 10 shows that 
the probability n for n = 2,4 configurations goes to zero 
at U > 700 meV, signaling the localization of electron on 
fullerene sites. Thus we obtain a critical ratio U/w = 1.4 
which is smaller than predi cted for assumed perfectly de¬ 
generate LUMO bandThis, however, does not imply 
automatically an insulating state for K3C60 since the up- 
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per recent estimate for U in this fulleride is 750 meVj^ 
and the actual val ue can be significantly lower as dis¬ 
cussed below (Sec. VIB). On the other hand, it would 
be incorrect to view the JT effect in the LUMO band as 
simply leading to its enlarging (Fig. Hit) which increases 
the critical U within (enlarged) single-band picture. As 
a matter of fact, the electron correlation and the metal- 
insulator transition in fullerides develops mainly in the 
middle adiabatic subband. The role of the middle band 
in the Mott-Hubbard transition can be qualitatively re¬ 
produced by single-band model. The band energy of the 
single-band model, which includes only one of the tiu or¬ 
bitals (the corresponding DOS is shown in Fig. [T^), is 
obtained as Et = —74.7 meV. Using the formula of the 
critical U for non-degenerate band within Gutzwiller’s 
approximation,l^we obtain U = 8|F^t| = 598 meV, which 
is close to U ^ 700 meV for AsCeo obtained in the 
present work (Fig. 10). The latter is larger than the es¬ 
timate for the single-band model by about 100 meV due 
to the remaining hybridization of the split bands. From 
this analysis, one may conclude that the Mott-Hubbard 
transition mainly develops in the middle band. 

An important issue is the accuracy of the calculated 
ground state energy. We used here a six parameter 
Gutzwiller ansatz (19) in combination with Gutzwiller 


approximation for the calculation of total energy (20). 
For comparison, Gunnarsson et alW^ used a single¬ 
parameter (conventional) Gutzwiller ansatz but calcu¬ 
lated the total energy without approximation within vari¬ 
ational Monte Garlo (VMG) approach. Gomparison with 
exact results obtained for small clusters of Geo via exact 
diagonalization has shown that VMG reproduce the ex¬ 
act total energy with accuracy of 0.1 % (see Table 7.1 
in Ref. l37|) . i.e., few meV of the total energy per one 
Geo in fullerides. The deviation from exact total energy 
will be certainly larger in the case of Gutzwiller approx¬ 
imation applied here, however, it will not cover the gain 
of the total energy due to JT splitting/orbital dispro¬ 
portionation amounting many tens of meV (Fig. [^). 
Thus the main conclusion concerning dynamical JT in¬ 
stability in fullerides seems to be unaffected by this ap¬ 
proximation. This is further corroborated by the fact 
that the Gutzwiller wave function used in the present 
work is more flexible in variational sense than in the con¬ 
ventional Gutzwiller ansatz. Indeed, even in the case 
of degenerate LUMO band, conforming to cubic symme¬ 
try, the ansatz (19) involves two projection parameters. 


These are An controlling the population of configura¬ 
tions UmAcr^mA-cr: ^ud An Controlling the population 
of configurations nmAcr^mA'cr', A 7^ A'. That these are 
the only independent parameters allowed by the cubic 
symmetry can be understood if one generalizes the form 
(19) to arbitrary LUMO basis on fullerene sites. Al¬ 


though, in general, UmAa in Eq. (19) are replaced by 


elements of one-particle density matrix, only the diago¬ 
nal remains nonzero because of the cubic symmetry. The 
terms in exponential of Eq. (19) then become of the 
form Aaaf3i3nmaai^mf3(j' ^ in which the “elasticity” tensors 


Aaai3i3 will be characterized by only two independent 
parameters in the case of cubic symmetry.^!! The sec¬ 
ond variational parameter in the Gutzwiller wave func¬ 
tion changes drastically the description of Mott-Hubbard 
transition in the cubic band. Thus, in a conventional 
single-parameter Gutzwiller ansatz within the Gutzwiller 
approximation the critical U/w = 4 in the case of three¬ 
fold orbital degeneracy of sites.l^ This critical ratio is re¬ 
duced to 2 in the case of Gutzwiller ansatz applied here 
(Fig. [1^) which is m uch c loser to values obtained by 
Monte Garlo treatment It would be of interest to use 
in future the Gutzwiller ansatz for static and dynamic JT 
effect on sites proposed here as a trial functions in vari¬ 
ational (VMG) and diffusion (projection) Monte Garlo 
(DMG) methodJ^ which give more accurate description 
of ground state energy. 


B. Parameters of the LUMO model 

Another important aspect concerns the values of rel¬ 
evant parameters for the model description of LUMO 
band in fullerides. Given the large number of such pa¬ 
rameters, their accurate knowledge is of primary impor¬ 
tance for realistic description of electronic properties of 
fullerides. Recently, it was proven that the DFT cal¬ 
culated vibronic constants of Ggg with a hybrid B3LYP 
functional compare well with those extracted from pho¬ 
toemission spectroscopy.!^ Then the vibronic constants 
for the GgQ anion and the exchange parameter calcu¬ 
lated within the same DFT functional should be reliable 
as well. Goncerning the transfer Hamiltonian, the pa¬ 
rameters of nearest-neighbor and next-nearest-neighbor 
tight-binding model (Table reproduce well the disper¬ 
sion of the LUMO bands calculated within LDA (gener¬ 
alized gradient approximation, GGA) (Fig. [l^ . It was 
shown by GW calculation that the interband electronic 
interaction can enhance the LUMO (tiu) bandwidth in 
fullerites by 30 9^^ while the intraband interaction re¬ 
duces the LUMO bandwidth in A3 Geo!^ However, in the 
latter case the GW approximation is ill-def ined d ue to 
strong correlation effects in the LUMO band.l^^EH 

The parameter assessed with less certainty in the 
model Hamiltonian 0-0 is the intra-fullerene electron 
repulsion U\\. Recent calculations of this parameter by 
constrained random-phase approximation (cRPA]P^with 
GGA band energies and wave functions give in the low- 
frequency limit U|| < 1 eV for the series of A3G60, the 
smallest being = 820 meV for K3G6o!^It is interest¬ 
ing to note that the estimated Uy in Ref. [31] gives for all 
fullerides values < 1 eV, while former estimation s made 
for fullerite (pure Geo crystal) give larger values. 

This is explained by the fact that in A3 Geo fullerides the 
LUMO Wannier orbitals occupy larger volume due to 
hybridization with alkali atoms.^^In these calculations a 
non-interacting polarization function was used that ex¬ 
cluded polarization processes within the LUMO band. 
Including the latter, i.e., considering the full screening 
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in the non-interacting metallic regimeSS further reduces 
U\\ by ca one order of magnitude in the low frequency 
limitA similar strong effect of metallic screening (aris¬ 
ing from the LUMO band) was predicte d also for a sim¬ 
plified treatment of metallic polarization.ESEI One should 
note that RPA is generally not expected to perform well 
in the limit of strong correlation, as well as the GW ap¬ 
proximation mentioned above. To check the screening 
capability of the correlated LUMO band, Koch et a/.SSI 
did VMC and DMC calculations of an induced charge 
arising in response to a test charge for a threefold degen¬ 
erate LUMO model of AsCeo- They found that RPA 
performs surprisingly well till U/w < 2, and even at 
the point of Mott-Hubbard transition (corresponding to 
U/w ^ 2.5 in their model), when the LUMO electrons 
become localized, the screening charge is reduced by not 
more than 40 % with respect to RPA screening charge 
calculated for non-correlated LUMO band (Fig. 2 in 
Ref. [48|) . At the same extent is expected to be re¬ 
duced the screening of the electron repulsion parame¬ 
ter by the intra-LUMO band interaction, which means 
that this screening is significant in the entire metallic 
phase of fullerides and should be taken into account for 
realistic assessment of U||. A rigorous model for the 
LUMO band, in which other degrees of freedom are ex¬ 
cluded, requires frequency-dependent electron repulsion 
parameters.!^ From it an effective static Hubbard model 
(involving a frequency-independent U) can be derived by 
fitting the self-energy in a low-frequency domain, not ex¬ 
ceeding the width of uncorrelated LUMO band.® Note 
that the derivations of the static U||, U± in Eq. Q 
should be done self-consistently with the derivation of 
the ground state following the iterative equations (j^- 

Experime ntally U can be assessed from Auger 
spectroscopy.SSSSl xhe estimates ar e 1.4 ± 0.2 eV for 
pure Ceo and band insulator and 0.6 ± 0.3 

eV for metallic KsCeo- The much smaller value of U 
in KsCeo reflects most probably the additional strong 
screening from half-filled LUMO band.® An instructive 
example of the sensitivity of U on intra LUMO-band 
screening is offered by non-cubic fullerides NHsKsCecP^ 
and (CH3NH2)K3C6 o!^ Contrary to the parent K3C60 
fulleride, which is a metal, these compounds are antifer¬ 
romagnetic Mott-Hubbard insulators. The main effect of 
spacers, NH3 and CH3NH2, respectively, is the removal 
of degeneracy of three LUMO orbitals on fullerene sites, 
which apparently reduces the orbital U/w from its value 
in cu bic K 3C60 predicted for threefold degenerate LUMO 
band^® causing the Mott-Hubbard transition. Manini 
et a/.® have checked this possibility via DMET calcu¬ 
lations of a model twofold degenerate band and found 
that the calculated splitting of the tiu orbitals is indeed 
sufficient to induce the Mott-Hubbard transition.® How¬ 
ever, the persistence of strong JT distortions in metallic 
fullerides, established in this work, calls for another inter¬ 
pretation. The crystal anisotropy induced by the spacers 
will enhance the splitting of the LUMO bands (Eig. Elf) 


thus reducing the intra-LUMO band screening of electron 
repulsion. This results in the increase of U\\ and U±_ in 
Eq. which is the reason why the non-cubic fullerides 
are Mott-Hubbard insulators. 

To conclude this part, several theoretical arguments 
and relevant experimental data argue in the favor of non- 
negligible intra-LUMO band screening of Hubbard U pa¬ 
rameter, which is thus expected to be well below 1 eV. 
The latter is also a necessary condition for metallicity of 
fullerides in the presence of strong JT distortions (Eig. 
10 ). 


Polaronic effects 


The simple form of dynamical vibronic wave func¬ 
tion (58) was derived under two simplifying assumptions. 
Eirst, the polaronic effect was neglected which seems to 
be justified for fullerides. Indeed, while the static JT 
energy of C^q and C^q is larger than in C^q by an 


^2 
-^60 

P(l) - 


amount /2 = 50.2 meV, the difference in 

the gain is compensated by the loss of stabilization en¬ 
ergy from the dynamical contribution. In the strong vi¬ 
bronic coupling limit, the dynamical JT contribution of 
Cgo is Zhjj/2 and that of Cgg Cgg is lEo because the 
trough is three-dimensional in the former case and two- 
dimensional in the latter. Then 1/3 of the dynamical 
contribution of Cgg is lost if JT relaxation accompanies 
the electron transfer. Using the data from the numerical 
diagonalizations, the loss of the dynamical JT contribu¬ 
tion is estimated as —F^djt/ 3 ~ 30 meV. Therefore, the 
binding energy of JT polaron (the energy gain arising 
from full JT relaxation) is AE = —Ej}^ — Edjt/ 3 ~ —20 
meV. Compared with the total JT stabilization energy of 
« 240 meV the JT polaronic effect appears to be small. 
One should take into account that the JT polaronic ef¬ 
fect is accompanied by the Eranck-Condon reduction of 
the band energy, which means that the JT polaron will 
only show up when the band energy is reduced by corre¬ 
lation effects under 20 meV, i.e., close to Mott-Hubbard 
transition. On the other hand the stabilization energy of 
one electron after total symmetric fullerene distortions 
does not exceed 20 meV, i.e., is negligible either.® In 
these estimations the relaxation due to displacements of 
alkali atoms has not been included, which is unimpor¬ 
tant for A3C60 but can be significant in insulating A4C60 
and AeCeo-!^ As for the second assumption of weak hy¬ 
bridization of the bands belonging to different tiu orbitals 
(Fig. §, it seems to be only justified in the strongly cor¬ 
related limit. When it is not the case, the 0-dependence 
of the coefficients U^ii in Eq. ( [4^ cannot be neglected 
and ultimately the rotations of JT deformation on differ¬ 
ent fullerene sites (Eq. ( [5^ ) cannot be separated. This 
means that in metallic CS3C60 the rotation of JT defor¬ 
mations occurs independently on different fullerene sites, 
while in K3C60 these rotations are more probable to be 
correlated. In the latter case the wave function (58) does 


not represent a close solution and should be rather con- 
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sidered as a variational function which nevertheless will 
correspond to lower total energy than the static JT solu¬ 
tion (ra. 


D. Summary 


0(Fig.i. Using the unit vectors of fee lattice, 
a a 

— 2 ^z) 1 ^2 — 2 ^x) "I 

as = I (©a, + ey ), (Al) 


The main achievements of this work can be summa¬ 
rized as follows: 


the displacements Am of the nearest neighbor sites from 
site m are written as 


1. We have developed an approach for the investiga¬ 
tion of correlated JT metals based on self-consistent 
Gutzwiller approximation. 

2. The concomitant treatment of JT effect and elec¬ 
tron correlation in metallic fullerides AsCeo proves 
the existence of dynamical JT instability in their 
ground state. The JT distortions arise due to 
strong reduction of the band energy by electron 
correlation effects and achieve an amplitude close 
to the value in a free Cgg ioii- 

3. The JT instability induces strong overall enlarge¬ 
ment of the uncorrelated LUMO band and its split¬ 
ting in three components corresponding to individ¬ 
ual adiabatic orbitals on fullerene sites. The results 
call for reconsideration of the role played by orbital 
degeneracy in the physics of metallic fullerides. 


Am = (a 3 ,ai - a 2 , -as, -(ai - a 2 ), 

3-1) ^2 — 3-35 —3i, —(a 2 — as), 

32,33-ai,-a 2 ,-(as-ai)). (A2) 

The next nearest neighbors are displaced by vectors 

Am = a(ea,,-e3;,ey,-ey,e^,-e^). (A3) 

Here, a is the lattice constant of a simple cubic lattice 
and correspondingly are unit vectors directed 

along tetragonal x^y^z axes (Fig. [^. The tight-binding 
Hamiltonian has the form: 

= E E + E E , (A4) 

m Act m cr 


4. JT distortions together with electron correlation 
induce disproportionation of electron density be¬ 
tween subbands corresponding to different adia¬ 
batic orbitals on fullerene sites. Besides the JT 
splitting there is also a bielectronic contribution to 
the separation of these subbands which vanishes in 
the limit of strong correlation. Importantly, the 
orbital disproportionation does not exist as a pure 
electronic low-symmetry instability in the absence 
of JT effect on fullerene sites {g = 0), in which 
case the correlated LUMO band will have a perfect 
cubic symmetry for any U. 

Finally, we note that a similar analysis can be applied to 
other correlated metals with JT active sites. 
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Appendix A: Tight-binding parametrization of the 
LUMO band structure of KaCeo 

We assume that all Geo’s in fee KsGeo lattice are 
equally orientated in a similar fashion shown in Figure 


where the nearest-neighbor part is 


= tl 
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FIG. 13. (color online) Band structure of fee KsCeo in eV. 
The red and blue points are obtained from DFT calculations 
and model Hamiltonian, respectively. The symmetric points 
{kx,ky,kz) are as follows: F = (0,0,0), X = (21^ja, 0,0), IF = 
(27r/a, 7r/a, 0), K — (37r/2a, 37r/2a, 0), L — (njjja)^ 
U = {27vla,7vl2a,7rl2a). 


TABLE I. LUMO level e (eV), transfer parameters, and band 
width w (meV) for KsCeo- 


e ti 

t2 

ts 

U 

tb 

te 

ty W 

5.066 43.3 

-31.9 

-6.2 

-16.6 

-9.6 

-2.0 

2.7 505.9 


and the next-nearest-neighbor part is 


-^mcr — ^5 {^m.-\-aGxXa^n 


F Cjn— 


m-aGxXa'-^mxa 


mya 


F Cjn— 


m—aeyya'^mya 


F CYn-\-aGzZa^mza F ^m—aezZa^mza^ 

F ^6 {^m-\-aGxycr^mya F Cm—ae^j^/crCmycr 
F CYn.-\-aeyZakniza F Cm —ae^^^crCn 


^mza 


Cm+ae;; 


xa'-^mxa 


F tj (Cm+ae^2;crCm2;cr Fc 
F Cm+ae^cccrCmxcr F Cm—ae^a^erC] 

F Cm+ae^ycrCmycr F Cm—ae^j^crCmycr) : 


m—ae^ a:cr'-ma:cr / 
■m—aGxZa^mza 
mxa 


(A6) 


Appendix B: Gutzwiller reduction factors for the 
LUMO bands and bielectronic energy in AsCeo 


To derive the form of Eq. (20), we apply the 


Gutzwiller’s approxi matio n extending the one for the 
nondegenerate band.^^^^ Within the Gutzwiller’s ap¬ 
proximation, physical quantities are described in terms 
of the probability z/p that one-site electron config uratio n 
F appears in the Gutzwiller’s wave function 

The occupation number Ux in spin orbital xa is de¬ 
scribed by the probabilities z/p as follows: 

2 ^ ‘^^XX F ^J^xy F 4z^ zx 


z 

F + 2Uzxz F ^T^zxx F ‘^T^xyy F 8l^xyz 
F ‘^l^xyxy F ‘^T^zxzx F SZ/'. 

+ 4z/. 


^xyzx F ^T^xyzy F ^T^xyzz 
^xyzxy F ‘^T^xyzyz F ^T^xyzzx F ‘^T^xyzxyz\-) (^t) 


where, 1/2 is due to the spin degrees of freedom, and 
A and A (A = x^y^z) indicate spin orbitals (A f) and 
(A I), respectively. Since we consider the metallic phase, 
z^p does not depend on the spin part of F. For example, 
^x = ^x- '^y and Uz are obtained by cyclic permutation 
of the indices (x^y^z) in Eq. (Bl). 

The Gutzwiller’s reduction factors Qxx appearing in Eq. 


(21) is given by 


Qxx 


1 


^x^xx + 2y ^y^x 


a yy'^yj'^x ' y'^x'^xx ' ^yj'^y'^xy 

-rix) 

F 2-\J'Uz^zx F ‘^yjT^yx^xyx F ‘sj^zz^zxz F zx^zxx 
F yj^yy^xyy F ^yj^yz^xyz F y/j^yxy^xyxy 
F \J^zzx^zxzx F ^ yzx^xyzx F 2 y^l^yzy^xyzy 

F ‘^yj^yzz^xyzz F ‘^yj^yzxy^xyzxy F -\J^yzyz^xyzyz 

F ‘^yjl^yzzx^xyzzx F s/l^yzxyz^xyzxyz^ •> 

where F = 0 means the configuration with no electron. 
Qyy and qzz are obtain ed b y cyclic permutation of the 
indices {x^y^z) in Eq. (B2). For gAA'(A 7^ A'), following 
relation holds: 


QXX' = V^AA^A'A'- (B3) 


The bielectronic energy is 


In Eq. (A5), Am^ indicates ith nearest neighbor. 

The DFT calculation of the band structure of KsCeo 
was performed using Quantum ESPRESSO 3.0 pack¬ 
age with the psendopotentials C.pbe-mt gipaw.UPF and 
K.pbe-mt fhi.UPF.I^ The lattice constant of KsCeo was 
taken from Ref. CT] and the structure of Ceo of Rof- M 
was used. 


The band structures from the DFT calculation (red) 
and the fitted tight-binding Hamiltonian (blue) are 
shown in Fig. The transfer parameters derived from 
the DFT calculation are tabulated in Table [ij The present 
values are close to the recent estimates with optimized 
structure.!^ 


-E'bi — ^11 {^XX F l^yy F l^zz) 

{U± — J}i/2){4:l/xy F 4:l/yz F ^l^'zx) 

F (3Ux F Fh)(2 ^xyx F 2lXyzy F ^J^ZXZ 
F ‘^J^zxx F ‘^J^xyy F 2Vyzz) 

F {3U± - 3J}i/2)8iXxyz 
F {QU±_ F 2Jh)( ^xyxy F ^yzyz F ^zxzx^ 

F {6U±_ — J}l/2){AUxyzx F ^T^xyzy F ^^xyzz^ 

F {}-^U±_)(2Uxyzxy F ‘^^xyzyz F ^^xyzzx) 

F _\_)Uxyzxyz' (B4) 


The Coulomb contribution (67) to the subband energy 
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level is given by 

_U 

^hi,x — [\^xx “1“ ‘^^xy ‘^^zxj 

'^x 

H- 2 i^Uxyx H“ ‘^^zxx ^xyy ^zxz H“ ^^xyz^ 

H“ 3 {yxyxy “1“ ^zxzx H“ ^^xyzx H“ ‘^^xyzy ‘^^xyzz) 

H- 4 i^Uxyzxy “1“ ‘^^xyzzx H“ ^xyzyz) 

“ 1 “ ^^xyzxyz] • (^ 3 ) 


Here, we choose the ordered JT distortion (|13|). ebi,^ 


and 


ebi,;2 are obtaine d by cyclic permutation oFthe indices 
(x^y^z) in Eq. (B5). Aebi in Eq. (68) is obtained as 

Aebi = ^hi,z ^hi,X' 
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